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Mathematical  Modelling  of  photovoltaic  (PY)  modules  is  important  for 
simulation  and  performance  analysis  of  PV  system.  Therefore,  an  accurate 
parameters  estimation  is  necessary.  Single-diode  and  two-diode  model  are 
widely  used  to  model  the  PV  system.  However,  it  required  to  determine 
several  parameters  such  as  series  and  shunt  resistances  that  not  provided  in 
datasheet.  The  main  goal  of  PV  modelling  technique  is  to  obtain  the  accurate 
parameters  to  ensure  the  I-V  characteristic  is  closed  to  the  manufacturer 
datasheet.  Previously,  the  maximum  power  error  of  calculated  and  datasheet 
value  are  considered  as  objective  to  be  minimized  for  both  models.  This 
paper  proposes  the  PV  parameter  estimation  model  based  minimizing  the 
total  error  of  open  circuit  voltage  (V0cX  short  circuit  current  (Isc)  and 
maximum  power  (PMax)  where  all  these  parameters  are  provided  by  the 
manufacturer.  The  performance  of  single-diode  and  two-diode  models  are 
tested  on  different  type  of  PV  modules  using  MATLAB.  It  found  that  the 
two-diode  model  obtained  accurate  parameters  with  smaller  error  compared 
to  single-diode  model.  However,  the  simulation  time  is  slightly  higher  than 
single-diode  model  due  extra  calculation  required. 
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1.  INTRODUCTION 

The  contribution  of  renewable  and  green  energy  sources  such  as  solar  photovoltaic  is  increased 
every  year  and  provide  clean  energy  to  community.  Therefore,  it  required  an  accurate  modelling  of  PV 
simulator  in  order  to  investigate  the  performance  of  PV  system  for  different  environment.  It  also  can  be  used 
to  evaluate  the  performance  of  converters  including  maximum  power  point  tracker  (MPPT)  as  well  as  system 
behavior  under  grid-connected  and  off-grid  [1],[2].  Basically,  PV  module  offers  solar  energy  conversion  into 
electrical  energy  by  using  solar  cell  [3].  The  Bell  Laboratories  revealed  that  silicon  had  photoelectric 

attributes  and  rapidly  developed  the  silicon  solar  cells  by  achieving  6  percent  of  its  efficiency  and  the  former 

satellites  were  the  elemental  use  for  these  first  solar  cells  [4] . 

A  PV  cell  is  a  semiconductor  diode  which  when  exposed  to  light  can  produce  charge  carriers  in 
order  to  generate  electricity.  If  the  PV  cell  is  connected  to  an  external  circuit,  current  will  flow  through  the 
circuit.  However,  a  single  PV  cell  produces  current  between  3  and  4  A  and  a  voltage  approximately  0.5  V.  It 
generates  very  low  power,  0.5  to  1.5  W  if  crystalline  silicon  cell  is  used  [5].  Hence,  series  or  parallel 
arrangements  of  cells  are  made  to  obtain  a  desired  electrical  voltage  and  power.  The  photovoltaic 

technologies  are  included  of  photovoltaic  cell,  module  and  array  devices  that  convert  light  energy  into 

electricity  directly.  PV  cell  is  made  up  of  semiconductor  materials  and  the  most  common  types  known 
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commercially  are  the  silicon,  mono -crystalline  and  poly-crystalline  [6].  The  number  of  the  PV  cells  are 
connected  in  series  to  increase  the  output  voltage  of  the  PV  module  [7].  While,  multiple  strings  of  PV 
modules  are  used  to  obtain  desired  power  the  practical  purposes  [8]. 

The  electric  circuit  based  on  the  PV  modeling  is  further  classified  into  three  types  of  model  such  as 
ideal  diode  model,  single-diode  model  and  two-diode  model.  The  simplest  model  is  the  ideal  diode  model 
which  comprises  of  a  single  diode  and  irradiation  dependent  current  source  [9] .  The  performance  of  the  ideal 
model  is  improved  by  adding  one  series  resistance  and  one  shunt  resistance,  it  is  generally  known  as  the 
single-diode  model  [10],[11].  Although  the  single-diode  model  is  more  popular  for  PV  modeling,  but  it  has 
weaknesses  such  as  it  neglects  recombination  loss  in  PV  cell  depletion  region  [12]. 

To  improve  the  accuracy  of  the  single-diode  model,  two-diode  model  has  been  introduced  by  adding 
another  diode  to  the  circuit  [13].  At  low  irradiance  and  temperature,  the  two-diode  model  offers  much  more 
accurate  curve  characteristics  as  compared  to  the  single-diode  model  [14].  The  single-diode  model  can  obtain 
fast  results  due  to  less  complex  equation  as  compared  two-diode  model.  But,  the  two-diode  model  provides 
more  accurate  characteristics  for  different  PV  modules  longer  iterations  [15]. 

In  PV  system,  Maximum  Power  Point  Tracker  (MPPT)  is  used  to  ensure  the  PV  module  operated  at 
the  Maximum  Power  Point  (MPP)  [16],  thus  total  power  output  from  PV  can  be  fully  utilized  by  the  system 
[17].  Various  MPPT  algorithms  have  been  proposed  in  order  to  increase  the  PV  system  performance  which 
includes  perturb  and  observe  (P&O),  open  circuit  voltage  (OCV),  short  circuit  current  (SCC)  and  Incremental 
Conductance  (IC)  [18], [19].  Therefore,  it  required  accurate  PV  simulator  to  test  MPPT  performances.  This 
highlighted  the  important  of  parameter  estimation  of  PV  module.  In  [13],  the  two -diode  model  is  used  and 
the  parameters  have  been  estimated  based  on  based  on  the  maximum  power  error  of  calculated  and  measured 
value  provided  by  manufacturer. 

In  this  paper,  an  improved  modelling  technique  for  the  single -diode  and  two-diode  models  are 
proposed  based  on  total  error  of  open  circuit  voltage,  short  circuit  current  and  maximum  power  (as  provided 
by  datasheet).  The  accurateness  of  the  model  has  been  validated  on  three  different  PV  modules  with  different 
I-V  characteristics. 


2.  PV  MODELLING 

This  section  explains  the  details  mathematical  modelling  of  single-diode  and  two-diode  models  for 
photovoltaic  modules.  The  determination  of  PV  parameters  for  modelling  and  simulation  according  to 
manufacturer  data  are  also  discussed. 

2.1.  Single-diode  Model 

Single-diode  model  is  a  simple  model  for  PV  module  consists  of  several  parameters  such  as  light 
generated  current  (/PV),  series  resistance  (. Rs ),  parallel  resistance  (. RP )  and  ideality  factor  (a).  The  equivalent 
circuit  can  be  modelled  as  in  Figure  1 . 


The  mathematical  formulation  of  current  and  voltage  relationship  can  be  written  as  follows  [10]: 
^  =  I  pv  ~  I d  ~  I sh  (1) 


I  Ipv  Iq 


exp 


rV  +  IRs  x 

v  aVT  j 


-1 


rV  +  IRs  x 

v  Rp  j 


(2) 
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where  IPV  is  the  current  generated  by  the  incidence  of  light,  I0  is  the  reverse  saturation  current,  VT  is  the 
thermal  voltage  of  the  diode,  Ns  is  the  cells  connected  in  series,  q  is  the  electron  charge  (1.60217646  x  10~19 
C),  k  is  the  Boltzmann  constant  (1.3806503  x  10‘23  J/K)  and  T  is  the  temperature  of  the  p-n  junction  in  K.  It 
is  important  to  mention  that  diode  ideality  factor,  a  is  a  constant  which  depends  on  PV  cell  manufacturer.  The 
equation  of  PV  current  as  a  function  of  the  temperature  and  irradiance  can  be  written  as  [10]: 

^v=7r-taw+*i AT)  (4) 

^ STC 

where  IPVstc  is  the  light  generated  current  at  STC  in  Ampere,  AT  is  the  T  -  TSTC ,  ( TSTC=25°C )  in  Kelvin,  G 
is  the  surface  irradiance  of  the  cell,  GStc  is  the  irradiance  at  STC  equal  to  1000  W/m2  and  K7  is  the  short 
circuit  current  coefficient  provided  by  the  manufacturer.  The  saturation  current  can  be  calculated  as 
follow  [10]: 


(lSCm+K,AT) 

exp 

>ocS7C  +KvAT~ 

-1 

1 

_ 1 

where  Kv  is  the  open  circuit  voltage  coefficient,  IscSPc  is  the  short  circuit  current  at  STC  in  Ampere  and 
Vocstc  is  the  open  circuit  voltage  at  STC.  These  parameters  provided  by  manufacturer  in  datasheet. 

2.2.  Two-diode  Model 

The  equivalent  circuit  for  two -diode  model  shows  in  Figure  2.  The  mathematical  equation  for  I-V 
relationship  can  be  derived  from  this  circuit  as  shown  in  (6)  and  (7)  [13]. 


vw 


V 


Figure  2.  Two-diode  model  circuit 
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exp 


ij  { V  +  IRs ) 
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where  I0i  and  I02  are  the  reverse  saturation  currents  of  diodes  1  and  2,  respectively.  VT1  and  VT2  are  the 
thermal  voltages  for  diodes  1  and  2,  respectively.  a7  and  a2  are  the  diode  ideality  factor  for  diodes  1  and  2, 
respectively.  The  I02  term  in  equation  (7),  compensate  the  recombination  loss  in  the  depletion  region  [13]. 
The  equation  for  PV  current  as  a  function  of  the  temperature  and  irradiance  as  in  (4)  [13]: 

I pv  (^pvstc  +^1^7")  (8) 

^ STC 
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The  equations  for  both  reverse  saturation  currents  IQi  and  I 02  are  set  to  be  equal  in  magnitude  and 
the  solution  can  be  achieved  analytically  using  the  following  equation  [13]: 


I  o\  1 02  I o 


(/„  +K.A7-) 


exp 


VocSK  +KvAT 


^  ax  +  a2  ^ 


(9) 


2.3.  Parameter  Estimation  of  PV  Module 

The  values  of  RP  and  Rs  are  achieved  through  iteration  by  calculating  them  simultaneously,  in  order 
to  match  the  calculated  peak  power  and  experimental  peak  power,  by  iteratively  increasing  Rs  value  while 
simultaneously  calculating  RP  value.  Under  the  maximum  power  condition,  the  equation  of  RP  can  be  derived 
from  (7)  and  expressed  as  [13]: 


RP  = 


(VMP  +  JMPRS  ) 


(/  _/  _/  \_'rMAX,e  / 

l/PV  1D\  1 D2  )  /V 

/  *  A/I 


The  initial  conditions  for  both  resistances  are  as  follows  [13]: 


Rso  —  0  ’  Rpo  — 


vu 


SC  ^MP  J 


V  -V 

v  OC  y  MP 


(10) 


(11) 


The  implementation  steps  for  determination  of  PV  module  parameters  based  on  single -diode  and 
two-diode  models  are  summarized  in  Figure  3.  The  manufacturer  datasheet  provided  the  important  points  of 
I-V  characteristic  such  as  ISc>  Voc  and  PMax  as  shown  in  Figure  4.  Therefore,  several  parameters  value  must 
be  estimated  in  order  to  model  the  PV  modules  based  on  the  equivalent  circuits  in  Figures  1  and  2. 

Previously,  the  estimated  parameters  of  PV  module  has  been  calculated  based  on  the  maximum 
power  ( PMax )  only  where  the  error  of  Voc  and  ISc  are  neglected  in  objective  function.  Therefore,  this  paper 
proposed  the  total  error  of  measured  and  calculated  values  of  ISc,  Voc  and  PMax  to  be  minimized  in  order  to 
obtained  the  accurate  PV  parameters.  The  total  error  is  formulated  as  follows: 

Error  =  |(^AX  c  —  Pmax ,e )  (Voc,c  ~ ^ OC,e )  SC,c  ~  ^ sc,e )|  (1^) 


where,  P max,c~P max,&  is  the  maximum  power  error  of  calculated  and  experimental  values  (datasheet),  Voc,c  ~ 
V oc,e  is  the  open  circuit  voltage  error  of  calculated  and  experimental  values  and  ISCfC  -  lSc,e  is  the  short  circuit 
current  error  of  calculated  and  experimental  values.  Newtown  Raphson  algorithm  has  been  used  to  minimize 
total  error  to  determine  the  best  value  of  RP  and  Rs  as  shown  in  Figure  3. 
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Figure  3.  Algorithm  for  minimizing  the  total  error  of  calculated  and  measured  data 


3.  RESULTS  AND  DISCUSSION 

The  effectiveness  of  the  single-diode  and  two-diode  models  for  parameters  estimation  have  been 
simulated  in  MATLAB  software.  The  modeling  methods  defined  in  this  paper  are  validated  by  measured 
parameters  of  chosen  PV  modules  provided  in  the  datasheets  at  Standard  Test  Conditions  (STC)  and 
compared  with  estimated  value.  All  the  experimental  data  are  taken  from  the  manufacturer’s  datasheet  of  the 
selected  PV  modules.  Three  different  PV  modules  have  been  considered  for  the  verification  purposes  such  as 
BP  Solar  MSX-60  [20],  Kyocera  KC200GT  [21]  and  Shell  S36  [22].  The  PV  module  specifications  are 
tabulated  in  Table  1. 
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Table  1.  Manufacturer  Data  of  Different  PV  Modules 


Parameter 

BP  Solar  MSX-60 

Kyocera  KC200GT 

Shell  S36 

Isc  (A) 

3.8 

8.21 

2.3 

Voc  (V) 

21.1 

32.9 

21.4 

Imp  (A) 

3.5 

7.61 

2.18 

VMp  (V) 

17.1 

26.3 

16.5 

Kv  (mV/°C) 

-80 

-123 

-76 

Ki  (mA/°C) 

3 

3.18 

1 

Ns 

36 

54 

36 

3.1.  Estimated  Parameters 

Based  on  the  data  provided  by  manufacturer  in  Table  1,  the  newton  raphson  algorithm  has  been 
used  to  minimize  the  total  error  in  (12).  The  main  objective  is  to  minimize  the  total  error  so  that  the  obtained 
parameters  of  PV  module  near  to  the  manufacturer  characteristic  as  in  Table  1.  Table  2  shows  the  total  error 
of  calculated  and  experiment  values  for  three  different  modules  at  STC.  It  shows  that  the  two-diode  model  is 
slightly  accurate  than  the  single-diode  model  where  the  total  error  is  smallest  for  three  type  of  PV  modules. 


Table  2.  Comparison  of  the  Total  Error  of  Calculated  and  Experiment  values 


Models 

MSX60 

KC200GT 

SHELL  S36 

Single-Diode 

Two-Diode 

Single-Diode 

Two-Diode 

Single-Diode 

Two-Diode 

Pmp  error 

0.1013 

0.1012 

0.0796 

0.0665 

0.0966 

0.0966 

Isc  error 

0.0096 

0.0044 

0.0190 

0.0062 

0.0036 

0.0013 

V0c  error 

0.1000 

0.1000 

0.1000 

0.1000 

0.1000 

0.1000 

Total  Error 

0.2109 

0.2056 

0.1986 

0.1727 

0.2002 

0.1979 

Simulation  time  (s) 

0.69 

1.66 

1.13 

1.76 

4.34 

5.11 

The  two-diode  model  can  give  significant  reduction  of  total  error  as  presented  in  Figure  5.  However, 
the  simulation  time  of  two -diode  model  is  slightly  higher  than  single -diode  model  due  to  additional 
computational  calculation  required  by  two -diode  model. 


MSX60  TypW&le  SHELLS36 


Figure  5.  Percentage  of  error  reduction  obtained  by  two -diode  model 


The  estimated  value  of  IQ ,  IPV ,  Rs  and  RP  for  the  single-diode  and  two-diode  models  are  presented  in 
Table  4  and  5.  The  I-V  curves  characteristic  of  the  two  models  for  S36  module  is  presented  in  Figures  6. 
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Table  4.  Parameters  for  the  Single-Diode  Model 


Parameter 

BP  Solar  MSX-60 

Kyocera  KC200GT 

Shell  S36 

a 

1.3 

1.3 

1.2 

Io(  A) 

9.09  x  10‘8 

9.82  x  10'8 

1.68  x  10'9 

Ipv  (A) 

3.80 

8.21 

2.30 

RP  (Q) 

164.54 

256.65 

1516.38 

Rs  (£0 

0.12 

0.20 

0.82 

Table  5.  Parameters  for  the  Two-Diode  Model 


Parameter 

BP  Solar  MSX-60 

Kyocera  KC200GT 

Shell  S36 

<27 

1 

1 

1 

a2 

1.2 

1.2 

1.2 

Ioi=Io2  (A) 

4.70  x  1010 

4.12  x  1010 

2.05  x  1010 

Ipv  (A) 

3.80 

8.21 

2.30 

RP  (Q) 

116.53 

128.40 

806.36 

rs  m 

0.27 

0.30 

0.89 

Table  6  presents  the  comparison  of  the  propose  total  error  approach  and  PMax  error  approach  [9]  for 
PV  modelling  under  standard  condition  (25 °C  and  1000  W/m2)  of  I-V  characteristic.  The  comparison  is  based 
on  the  two-diode  model  in  order  to  highlight  the  effectiveness  of  proposed  approach.  It  found  that  total  error 
approach  can  provide  lower  PMax  error  approach  error  as  compared  PMax  error  approach. 


Table  6.  Comparison  of  PMax  Error 


Pmax  from  datasheet  (W) 


Method 


Pmax  Error  [19] 


Proposed  Total  Error 


Pmax=  36 


Estimated  Pm4x=35.9000  W 
Error=0.1000 


Estimated  PMax= 35.9034  W 
Error=0.0966 


4.  CONCLUSION 

The  modelling  of  PV  equivalent  circuit  based  on  single-diode  and  two-diode  models  have  been 
compared  in  this  paper  for  three  different  modules.  To  obtained  the  accurate  parameters  estimation,  the  total 
error  of  open  circuit  voltage,  short  circuit  current  and  maximum  power  of  calculated  and 
experiment/manufacturer  data  are  used  as  objective  function  to  be  minimized  using  Newtown  Raphson 
approach.  From  the  simulation  results  using  MATLAB,  it  shows  that  the  two -diode  model  is  better  than  the 
single-diode  model  in  term  of  total  error  minimization  where  it  can  be  reduced  up  to  13.04%  for  KC200GT 
module.  The  estimated  PMAX  produced  by  proposed  approach  is  significantly  lower  than  existing  approach 
where  all  the  i-V  characteristic  under  standard  condition  can  be  satisfied  efficiently.  On  the  other  hand,  the 
simulation  time  of  two-diode  model  is  slightly  higher  than  single-diode  model  because  it  required  extra 
calculation.  It  can  be  concluded  that  parameter  estimation  for  two -diode  model  based  on  total  error  reduction 
of  manufacturer  data  give  smallest  error  compared  to  single -diode  model. 
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